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Abstract 

We propose a new parametrization of the four-point vertex function in 
the one-loop one-particle irreducible renormalization group (RG) scheme for 
fermions. It is based on a decomposition of the effective two-fermion interac- 
tion into fermion bilinears that interact via exchange bosons. The numerical 
computation of the RG flow of the boson propagators reproduces the leading 
weak coupling instabilities of the two-dimensional Hubbard model at Van Hove 
filling, as they were obtained by a temperature RG flow in [l^] . Instead of 
regularizing with temperature, we here use a soft frequency fi-regularization 
that likewise does not artificially suppress ferromagnetism. Besides being more 
efficient than previous ./V-patch schemes, this parametrization also reduces the 
ambiguities in introducing boson fields. 

1 Introduction 

In the past decade fermionic Wilsonian renormalization group (RG) methods have 
been very successful in classifying weak coupling instabilities of the two-dimensional 
(t, t')-Hubbard model 0, 0, El, HE ESQ. The starting point in these methods 



is an exact functional equation. Besides being very useful for proving mathematical 
statements, this equation can be approximated to allow a direct numerical calcu- 
lation. In the case of weak coupling an expansion in the fermionic fields can be 
truncated by setting the six-point function to zero. For an appropriate Fermi sur- 
face geometry this one-loop approximation is also justified for intermediate couplings 
and scales [25]. 

Successes of the fermionic one-loop RG method include the explanation of the in- 
terplay of antiferromagnetic and superconducting instabilities in the Hubbard model 
for small next to nearest neighbor hopping — t 1 . Results were first obtained using 
Polchinskis version of the RG (32L I33I , the Wick ordered scheme @, 0], and the one- 
particle irreducible (1PI) scheme [131. 1 151] . Near half filling, where the Fermi surface 
is nested, antiferromagnetism was found to be the leading instability. For smaller 
fillings on the hole doped side, where the Fermi surface is more strongly curved and 
regular, d x 2_ y 2-wave superconductivity is dominating (although not present initially 
in the repulsive Hubbard model, such a term is generated as an effective interac- 
tion in the RG flow; specifically, it is induced by antiferromagnetic correlations). 
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In between the two regions the antiferromagnetic and superconducting tendencies 
mutually reinforce each other 0, This saddle point regime, where the main 
contribution comes from the saddle or Van Hove points, is interpreted as the weak 
coupling analogue of the Mott state, but not yet fully understood. 

In all works mentioned in the last paragraph a Fermi surface momentum or 
frequency cut-off was used. Although being conceptionally clear, this regulariza- 
tion suppresses small-momentum particle-hole fluctuations, as pointed out in [l2T |. 
Overcoming this drawback by using temperature as a flow parameter, the leading 
instability for larger hopping —t' £ [0.34i, 0.5t] was found to be ferromagnetism at 
Van Hove filling and triplet superconductivity away from Van Hove filling [ll|) E3] • 
This result was qualitatively confirmed by the same method in [18], by an inter- 
action RG flow [lol |. and a two-particle self-consistent Monte Carlo approach [sj]. 
The mutual suppression of ferromagnetic and d-wave superconducting tendencies 
decreases the critical scale at Van Hove filling for intermediate —t' by several orders 
of magnitude. Within numerical accuracy, the instability analysis of [ll|, [13] even 
suggests the existence of a quantum critical point between the ferromagnet and the 
superconductor. 

The one-loop 1PI RG flow used here is described by the evolution of an ef- 
fective four-point vertex function and a selfenergy. Even if the selfenergy is not 
taken into account, the flow equation for the four-point function is a non-linear 
integro-differential equation (see Section [2|). This equation has to be solved for a 
function with three independent momenta and three independent frequencies. For 
a numerical computation, as applied by the above mentioned studies, further ap- 
proximations are needed. By Taylor expansion and power counting it can be argued 
that the main contribution to low energy excitations comes from frequency zero. All 
above mentioned one-loop studies neglect the frequency dependence of the vertex 
function and evaluate the right hand side of the flow equation at frequency zero. 
In a next step the momentum dependence is discretized. Again by a low energy 
argument, the momenta of the vertex function are projected onto the Fermi surface 
in most one-loop studies. In a so-called TV-patch scheme (3, Is3] the Fermi surface 
is divided into N patches^] and the angular dependence of the vertex on each mo- 
mentum is approximated by a constant in each patch. Evaluating integrals as sums 
over patches, this leads to a coupled system of ~ iV 3 ordinary, non-linear differential 
equations, which can be solved numerically. In [l(| the momentum dependence is 
not projected onto the Fermi surface, so that full momentum space (— vr,7r] 2 has to 
be covered by patches. 

Despite the successes of the one-loop A^-patch schemes, the method shows room 
for improvement. - First of all, the selfenergy is neglected in most one—loop RG 
studies. Quasi particle scattering rates of the Landau Fermi liquid are calculated in 
BEE HE El, which gives a first hint of selfenergy effects. Deformations of the 



1 This only applies for N > 2. The extreme case are the two-patch models, where, contrary to 
TV-patch schemes, full momentum space is not taken into account. At Van Hove filling they expect 
the main contribution to come from two small patches around the Van Hove points [3. HR . 
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Fermi surface due to the real part of the selfenergy are taken into account in [15| . A 



new dynamical adjustment of the propagator [24( can help to include the selfenergy 
entirely. - A second problem is the observed flow to strong coupling, which limits 
the validity of the one-loop RG flow. The flow has to be stopped before all scales are 
integrated out. That is, the full model is not recovered in the flow since the infrared 
regularization is not removed completely. This drawback may be overcome with a 
modified one-loop scheme by allowing the flow to continue in a symmetry broken 
phase 17, 2f|. Using this modified scheme the exact solution of a mean-field model 



for superconductivity is obtained in the thermodynamic limit Furthermore, 
the scheme has been recently applied to the attractive Hubbard model [g]. Due 
to the emergence of an effective superconducting gap, which regularizes infrared 
divergencies, all scales can be integrated out. However, since the vertex function is 
not charge invariant anymore, the flow equations become very complicated. Besides 
the purely fermionic description of the RG flow, a bosonic continuation with mean- 
field theory [22| or partially bosonized RG flows [3, 2(3, 23,113] are promising methods 
to treat the symmetry-broken phases. 

We also address a third problem of the one-loop iV-patch schemes, which is 
the parametrization of the effective vertex function. Although the initial four-point 
vertex function of the Hubbard model is a constant in momentum space, a non- 
trivial and physically important momentum structure evolves in the one-loop RG 
flow. The approximations entailed by the projections in iV-patch schemes are not 
yet fully understood. A severe restriction on the possible choice of N is the high 
computational cost of solving ~ iV 3 ordinary differential equations, whose coeffi- 
cients have to be calculated as two-dimensional integrals in every integration step. 
For the same reason, including the frequency dependence is very costly. A more effi- 
cient parametrization of the four-point vertex function should identify the relevant 
processes. Their separation from irrelevant remainders can be guided by previous 
one-loop studies and the observation that the identification of the leading instabili- 
ties is determined by the asymptotic singular structure of the flow equation. 

In Section O we propose a new parametrization of the vertex function in the 
one-particle irreducible RG scheme. Guided by the singular momentum structure 
of the right hand side of the flow equation, we identify three channels with distinct 
singular momentum structures El All graphs of the flow equation are uniquely as- 
signed to a channel and the channels are general two-fermion interactions, so this 
is no approximation. However, since we define the channels based on their singular 
momentum structure, there are still ambiguities involved. This condition allows dif- 
ferent definitions, for example, how a constant term is distributed over the channels. 
Therefore the initial on-site interaction is not assigned to a channel but kept fixed 
in the flow explicitly. Nevertheless, the definition of the channels in terms of graphs 
is natural from the structure of the flow equation. 

In the thus defined channels the effective two-fermion interaction is expanded 



2 In a similar approach the effective vertex function of the single impurity Anderson model was 
recently written as a sum of three functions, each dependent on only one singular frequency [lr3 |. 
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in terms corresponding to fermion bilinears that interact via exchange bosons. For 
this study, we choose scale- and frequency-independent form factors for the fermion 
bilinears. In Section 01 we show that only a small number of terms is needed to cap- 
ture the essential features of the one-loop RG flow in the case where one channel is 
dominant. In the case of competing channels, we choose the parametrization guided 
by the results from A^-patch studies. Inserting this decomposition into the RG equa- 
tion and projecting the right hand side onto the coefficients of the expansions, the 
flow equations for the boson propagators are derived. In Section [6] we state these 
flow equations for general expansions with a finite number of form factors when the 
tails of the expansions are neglected. 

We test the proposed parametrization of the vertex function in an application 
to the (t, t')-Hubbard model at Van Hove filling and temperature zero. Because 
at this filling, the interplay between superconductivity and ferromagnetism is im- 
portant, it is essential to choose a scale-dependent regularization that does not 
artificially suppress small momentum particle-hole fluctuations [12 ] . Here we do 
not use the temperature flow, but instead choose a very mild infrared regulator that 
keeps these fluctuations, but provides enough regularization to have a well-defined 
flow equation, see Section [5j Further approximations are applied for the numerical 
implementation in Section [7J where we neglect the frequency dependence of the bo- 
son propagators and discretize their momentum dependence using step functions. 
Since this momentum dependence is mainly determined by the one-loop bubbles, 
there is good guidance in that. Also, we choose a small number of form factors. 
The numerical results are stated and discussed in Section [5] We obtain the same 



ordering tendencies of the Hubbard model at Van Hove filling as in [12] . 

In the remaining part of this introduction we describe the model and our no- 
tations. Consider the two-dimensional quadratic lattice T = Z 2 /LZ 2 with unit 
spacing, length L 6 2N, and periodic boundary conditions. Then the dual lattice 
or momentum space is given by the torus T* = ^Z 2 / '(27rZ 2 ) . Let a P)Cr and a PiCT be 
the fermion operators on the associated Fock space with momenta p S T* and spin 
a € {+, — } in quantization direction. The (t, i')-Hubbard model on T is defined by 
the Hamiltonian 

H[a\a}= e (p) a P ,^°W + U ^ a P i,+ a P2 - a P3 ,- a Pi+P2-p 3 ,+ (!) 

per* pi...p 3 er* 

with an on-site repulsion U > between electrons. The kinetic term of ([1]) emerges 
in a tight binding approximation. That is, the transfer integral, which describes 
electron hopping between different lattice sites in real space T, is set to t for nearest 
neighbors and —t' for next to nearest neighbors. In momentum space this corre- 
sponds to the dispersion relation 

e(p) = — 2t(cosp x + cos p y ) — At' cosp x cosp y — fj, , (2) 

where the chemical potential [i is already included in the definition of e(p) and TL, 
which is convenient for setting up the grand canonical partition function. 
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In the following t is set to one, that is, all quantities are measured in energy 
units of t. The next to nearest neighbor hopping is considered in the parameter 
range —t' E (0, |). It determines the shape of the non-interacting Fermi surface 
together with the chemical potential. For fi = At' the system is at Van Hove filling 
where the Fermi surface intersects the Van Hove points, where the gradient of the 
dispersion ([2]) vanishes. This causes a logarithmic divergence in the density of states. 



2 The Fermionic 1PI One-Loop RG 

We consider the one-particle irreducible (1PI) fermionic renormalization group (RG) 
flow in the one-loop approximation [25(. Since we restrict the RG flow to the 
symmetric phase, the effective interaction, which depends on the RG scale A, is 
assumed to be SU{2)— and [/(l)-symmetric, as is the initial Hubbard interaction. 
An effective two-fermion interaction with those symmetries can be generally written 
in terms of a four-point vertex function V in the following way 

Va.[*] = ^ J dp 1 ...dpiS(p 1 +p 2 -p3-Pi)VA(pi,P2,P3) (3) 

x ^2 Vv(pi)Vv(p2)Vv'(p3)tMp4) • 

6{ + ,-} 

Due to momentum conservation, the four-point vertex function Va is a function 
of three independent momenta and frequencies. Here we denote p = (po,p) with 
Matsubra frequency po and spatial momentum p. The integral J dp is shorthand 
notation for jj ^2 Po ^ja ^Jfi i n the thermodynamic limit L — > oo and at in- 
verse temperature (3. The initial on-site interaction of the (t, £')-Hubbard model is 
characterized by the constant vertex function Vh.m.(p1)P2>P3) = U. 

The flow equation for the effective vertex function in the 1PI RG scheme reads 

Va(pi, ■ • -Pa) = % P (Pi, ■■■P3) + 1^(pi, ■■■Ps)+ T^iPi, ■■■Ps), (4) 

where the dot denotes the derivative with respect to the scale A. The particle- 
particle contribution and the crossed and direct particle-hole contributions are re- 
spectively given by 

T PP (pi,---P3) = ~ J dp [-j K G{p)G(p 1 +P2~p)] V a {p 1 ,P2,p)Va(pi +P2~P,P,P3) 

2£(Pi,---flO = - j dp [m G (P) G (P + P3-Pi)] Va(pi,P + P3-Pi,P 3 ) 

x V A (p,P2,P + P3 - Pi) 
T p i(Pi,---P3) = J dp [j K G{p)G{p + p 2 -Pz)] 2VaOi,P + £>2 - P3,p) V A (p,P2,P3) 

~ V A (pi,p+P2 ~P3,Pl +P2 ~P3) Va(P,P2,P3) 
~ V A (pi,p + P2 -P3,P) Va(P,P2,P + P2 -P3) 
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The scale derivatives in the integrands act only in the square brackets. Since we 
neglect the selfenergy, the propagator G(p) = [ipo + e(p)]~ 1 XA(p) is the free propa- 
gator multiplied with an appropriate cutoff function. That is, although not denoted 
explicitely, it is understood throughout that the propagator G(p) depends on scale 
A. This flow equation is derived in [25| . 

Since the initial vertex function U of the Hubbard model obeys the following 
symmetries, the flow equation implies that the effective vertex function Va satisfies 

Va(pi,P2,P3) = Va(P3,P1 +P2 -P3,Pi) (PHS) _ 
Va(pi,P2,P3,) = V A (P2,P1,P1+P2-P3) (RAS). ° 

The first symmetry expresses invariance of the partition function as a functional 
integral under a change of integration variables tp i— > tip and tp i— > itp and is called 
particle hole symmetry (PHS). The second symmetry is directly inherited from ([3]) 
and accounts for remnants of the antisymmetry (RAS) of the Grassmann variables. 



3 Decomposition of the Effective Interaction 

The parametrization of the vertex function developed in this section is based on 
the observation that mainly the singular momentum structure of the RG equation 
determines the qualitative instabilities of the flow. That is, the parametrization of 
V should simplify the momentum dependence but should keep track of all possible 
singular contributions. If the vertex function is regular, the only momentum de- 
pendence that can change the singular behavior of the RG equation is the transfer 
momentum that propagates through the scale-derivative of the particle-particle or 
particle-hole bubble 

*PP/p h (0 = ^/ dpG(p)G(l T p). (6) 

In order to absorb this momentum dependence we introduce three additional chan- 
nels corresponding to three different transfer momenta in the RG equation respec- 
tively 

v A [*] = v h .m.[*] + v s A c [*] + y£[¥] + v£[*] ■ (7) 

The initial on-site interaction of the Hubbard model V h .m. [^] is kept in the parametriza- 
tion since there are ambiguities in assigning it to the other channels. It will remain 
constant in the flow and can be seen as a driving force. This does not mean that 
the on-site term contained in the full interaction ([7|) remains independent of scale. 
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The corrections to it are absorbed as contributions to the additional channels 



3 

V s A c [f] =-- ! dqdq'dl ^ c (q,q',l) £ (^cx^^/ - q)) - q' 

J J=0 

3 

Km = -\ f d^'dz ^(g, 9 ^oE(*(g)^ ) *(g+^))(w> a) *(g , -o 



Km = -\ j dqdq'dl ^(g,g^0(*(g)*(g + 0)(*(g / )*(g , " 0) , (8) 

where o - ^ are the Pauli matrices for J = 1,2,3 and <r(°) is the two-dimensional unit 
matrix. Each of the channels ([8]) is a general U(l) and SU(2) symmetric two-fermion 
interaction. We call them the superconducting, magnetic, and forward scattering 
channel respectively. These names are legitimate if the functions <3?sc> and ^ K 
are regular in their first two momentum indices and possibly singular only in the 
third (bosonic) momentum index I. Note that the $'s are zero at the beginning and 
are generated in the flow. In order to define their evolution we use the standard 
Fierz identity | J2 3 j=o ^^2^30-4 = S aia4 5 a20 - 3 to restate flHJ) in terms of a vertex 
function 

VA(pi,P2,P3) = U - $£ c (pi,P3,Pl +P2) + $m(j>1,P2,P3 -Pi) (9) 
I 1 
+ 2^mOi>P2,P2 -P3) - 2 $ k(P1 ; P2,P2 -P3) 

and insert it into the RG equation. The strong momentum dependence of the 
superconducting channel should be the sum of the two incoming particles p\ + 
P2 which is the transfer momentum of the particle-particle graph. Likewise the 
combinations p% — p\ and P2 — P3 are the transfer momenta of the crossed and direct 
particle-hole graphs respectively. So we define 

$sc(Pl>P3,Pi +P2) = -% P {PUP1,P3) 

®i(Pi,P2,P3-Pi) = T;i{p llP2 ,p3) (10) 

r p d h (pi,P2,P3) + T p Z( 



$k(J?1,P2,P2 -P3) = -27"h(pi,P2,P3) +^h(Pl,P2,Pl + P2 ~ P3) 



This decomposition of the interaction vertex into three channels is exact at one- 
loop level and the symmetries (PHS) and (RAS) are satisfied by each channel, that 
is 

^ c (q, q', I) = $£ c (l -q,l- q' , I) = $£ c (</ , q, I) (11) 

v', = ®wAq',q, -0 = ®wAq + i,q'- i, -0 • 

The first identity in each line corresponds to (RAS) and can be directly seen from 
dHj), whereas the second identity corresponds to (PHS). Both symmetries follow from 
the defining flow equations (fTUj) , Therefore ([9]) satisfies (RAS) and (PHS). 
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In summary, the definition of the channels is chosen such that each channel carries 
one critical momentum dependence. If the vertex function is still regular, all possible 
singular momentum dependence of the flow equation is absorbed by the channels. 
Furthermore, the decomposition satisfies (RAS) and (PHS). The assignment of the 
graphs, however, is not unique for special momenta. For example, a constant term 
can be freely distributed among the channels. Nevertheless, we have uniquely defined 
the channels assuming that the different classes of graphs with their corresponding 
transfer momenta form a natural decomposition. This is the only definition such 
that each channel absorbs one singular momentum and such that (if the transfer 
momentum becomes singular) the channels correspond to interactions of Cooper 
pairs, spin operators, and density operators respectively as written in (jHJ). 

Next we further decompose the interaction vertex by writing each channel as an 
interaction of two fermion bilinears interacting via an exchange boson. Naturally, 
the critical transfer momentum, which as a sum of two fermion momenta is a boson 
momentum, is chosen to be the momentum of the boson propagator. Due to (PHS), 
3>sc(<7> q'> = ^sc^'i OS an d we can expand 

*so(?.9',0= E D m n(l)f m ^-^fn^-^)+R sc (q,q',l), (12) 

such that Vsc[^] describes the interaction of Cooper pairs via the boson propagator 
D{1) up to a remainder term -R S c- Actually, $ sc as a function of its first two 
indices can be regarded as the kernel of a Fredholm operator. But since we are not 
able to determine the scale-dependent eigenfunctions, we expand in a given (scale- 
independent) finite set of form factors (/ n )neXsc anc ^ l eave a remainder. The form 
factors are orthonormalized on full momentum space (— vr,7r] 2 (divided by 4-7T 2 ) and 
depend only on the relative momentum of the Cooper pairs, representing different 
singlet and triplet lattice symmetries. In principle an expansion in momentum and 
frequency degrees of freedom is possible. For simplicity we take an approach where 
the form factors (in contrast to the boson propagator D) do not depend on frequency. 
This leaves an ambivalence in the definition of D such that the expansion (|12p is 
satisfied. We set the frequency corresponding to q and g' to y, which is a fermion 
frequency since Iq is a boson frequency. The leading instability will occur at Iq = 
so in that case qo = q' = which will give the main contribution. For (3 < oo this 
is an overestimate since the lowest possible (absolute) value of a fermion frequency 
is ±§. Alternatively for finite temperature, qo and q' could be set to ^ ± 5 and 

^ =F § respectively, with symmetrization over the sign. That is, for m,n E I ac , the 
coefficient D mn (l) in (fT2"j) is given by 

D mn (l) = J ^ ^ Uk-^Uk-^soiq^^l)^,^ • (13) 

With this choice of definition, (RAS) implies D mn (l) = D nm (l) = a m a n D mn (l) 
where f n (p) = cr n f n (—p) with a n either + or — . Therefore, D is a block matrix 
with separate blocks for singlet and triplet symmetry. 
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Likewise the magnetic and forward scattering channels are expanded. Now by 
(RAS), ^ m /k{q,q'J) = ^m/k(q',Q,~1), so we can write 

<f>i( q , q ',i)= ^W/m(q + ^)/n(q'4) + fi M(«,9'0 (14) 

m,neJMK 

$k(q,q',1)= E ^(0/m(q+|)/ n (q' -l)+R K (q,q',l) 

m,n£l MK 

with symmetric boson propagators M mn (l) = M nm (—l) and K mn (l) = K nm (—l) and 
remainder terms R M and R K . The functions (f n )n£i MK are again scale-independent, 
orthonormalized on (— ir, ir] 2 , and frequency independent. Although different expan- 
sions for the magnetic and forward scattering channel are possible, we choose the 
same for notational simplicity. The index sets I sc and I MK need not be the same. 
In analogy to the superconducting boson propagator we set 

M mn (l) = J ^2-^2 /m(q+ " D $ m(</, O^^k ( 15 ) 

/d 2 q d 2 q' 
(2tt) 2 (2tt) 2 ^ m ( q + s)/*^' ~ ^) $ K(g,g',0| 9o= _ii, ig , = ^ • 

With the boson propagators defined in this way the three expansions and the remain- 
der terms fulfill (RAS) and (PHS) separately. The magnetic interaction describes 
interacting spin operators. For example, a constant /s(p) = 1 describes a local 
spin operator and a further expansion in the s-wave channel generalizes to nearest 
neighbors and next to nearest neighbors and so on. A d-wave form factor in the 
forward scattering expansion describes a possible Pomeranchuk instability [a]. 

So far we have decomposed the interaction vertex such that boson propagators 
carry the critical momentum dependence of the right hand side of the RG equation 
at least as long as the effective vertex function is regular. If it can be shown that 
the remainder terms -R sc , R M , and R K remain regular or at least less singular than 
the boson propagators (even if the effective vertex function develops a singularity 
in the flow), then several benefits are gained for the analysis of competing insta- 
bilities. First, the decomposition of the interaction allows to identify qualitatively 
which instabilities are favored. After the flow is stopped, the dominant terms in 
the interaction can be decoupled by Hubbard-Stratonovich transformations, so that 
one can proceed with a bosonic flow. Furthermore, the parametrization of the flow 
is simplified, since no function of three fermion momenta and frequencies has to be 
studied, but only several functions of one fermion momentum and frequency. The 
boson fields involved have point singularities, which pose substantially less numerical 
effort than the extended Fermi surface singularities. 

4 The Superconducting Channel 

In this section we consider an example to illustrate our method. Let the Fermi 
surface be curved and regular. By this we mean that the Fermi surface does not 
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meet the Van Hove points and that Umklapp scattering is irrelevant. Then the 
particle-hole bubble is negligible compared to the particle-particle bubble, which 
develops a strong peak for small momentum and frequency. 

In a first step the particle-hole graphs are altogether neglected. The RG equation 
(jU) without selfenergy effects is then solved by the solution of a self-consistency 
equation 

V A (q, l-q,q') = V Ao (q, l-q,q')-f dp G(p)G(l - p)V Ao (q, I - q, p)V A (I - p, p, q') , 

(16) 

where G{p) = [ipo — e(p)]~ 1 XA(p) is the free fermion propagator multiplied with 
a regulator function xa(p)- In this section we choose a strict Fermi surface mo- 
mentum cut-off, that is, modes with |e(p)j < A are suppressed completely. If the 
vertex function is expanded in a complete set of A independent functions f n but A 
dependent coefficients Vj^ n (l) 

V A {q, I - q, q) = V A {q , I — q, q) = ^ v^ n {l)f m (± - q)f n {^ - q) , 



then the self-consistency equation (j!6f) is solved by v = A~ , where the scale- 
dependent matrix A mn (l) = (v Ao )mn(l) + ^^ n (0 is given by the initial condition 
and the particle-particle bubble 

$r(0 = / dpG(j>)G(l-p)Ul-p)f n (l-p). (17) 

We concentrate on the singular case 1 = 0. If the functions f n could be chosen such 
that 3>^ n (0) and (^ Ao )m},(0) are both diagonal, then the matrix ^4(0) could easily be 
inverted to give f„ n (0) = ^°(0)[1 + ^ (0)$^ p n (0)]- 1 5 mn . If there is an attractive 
channel in the initial condition, that is < for some n, the flow in this channel 
diverges for large enough f3 due to a zero in the denominator. On the other hand, 
the flow of is asymptotically free for all n with > (repelling). 

However, $™ ri (0) is not diagonal in general. Suppose that the functions f n do not 



depend on frequency. Following 28f] we change variables p = tt(E, 6) with E = e(p) 
and an angle 9 to obtain 

/ A ° tanh — r 
dE J de j(E,e)f m (-K(E,e))f n (7r(E,e)) 

A 

with Jacobian J(E,9). The ^-integral is a smooth function of E for a curved and 
regular Fermi surface. For Ao small enough the zeroth order of a Taylor expansion 
in E, that is a projection onto the Fermi surface, gives diagonality if the f n 's are 
chosen as Fermi surface harmonics [H. Therefore the main contribution of $™ n (0) 
can be made diagonal for a curved and regular Fermi surface. 
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This gives a good understanding of the flow in the superconducting channel (|12|) . 
In a diagonal expansion of dominant processes only form factors that give rise to 
positive boson propagators have to be taken into account. All other terms in the 
expansion are suppressed to zero and can be neglected since they will not influence 
the qualitative behavior of the flow. In practice one suspects that only the biggest 
positive coefficient plays a role in the flow to strong coupling (depending on the size 
of the corresponding bubble). 

For a non-diagonal expansion one cannot separate the irrelevant modes exactly. 
There will be a flow to strong coupling if det A = (for any finite expansion) at 
some scale. However, since diagonality holds approximately, there is a good chance 
to capture the singular behavior of A in a small matrix of well-chosen form factors. 

In a second step we take into account the full one-loop RG equation with all 
particle-particle and particle-hole graphs and start with the initial repulsive Hub- 
bard interaction. Due to the curved and regular Fermi surface and negligible Umk- 
lapp scattering assumed in this section, we only consider the superconducting chan- 
nel and decompose 

V A (q,l-q,q') = - £ D mn {l)f m {\ - q)/„(I - q') + U + R(q,l - q,q') . (18) 

m,n£lsc 

As discussed above we select only a few form factors that describe the supercon- 
ducting channel correctly and drop the remainder term i? sc of the expansion in the 
superconducting channel (|12p Fl The function R in (|18p arises from particle-hole 
graphs. We will study its influence on the flow in the superconducting channel and 
show that the Kohn-Luttinger effect is present in our method. That is, although 
starting at Ao with a repelling V\ (q,l — q,q') = U > 0, particle-hole terms will 
create an attractive superconducting interaction. 

Inserting the decomposition (fT8|) into the RG equation (JH) and projecting the 
particle-particle graphs according to (Tl3]) gives the flow equation for the supercon- 
ducting boson propagators for n,m G I sc 



Dmn(l) =cr n dfi(p,l-p) 



^ D ma( l )fa(^ ~ P) ~ S m , S U - a m (p,l) 

aeZsc 



(19) 



XI Aia'(0/a'(l - P) -8n,sU - On(p,l) 
a'eX sc 

with initial values D mn (l) = at scale Ao and with the contribution from the 
particle-hole channels a m (p,l) = f p^/m(| - <l)R(qJ ~ q,P)\ qo =i /2- We denote 
dn(p,k) = gjj- (G(p)G(k)) dp for the bubble integration and f m (— p) = cr m f(p) 
distinguishes between singlet and triplet symmetry. Further we assume that I sc 
contains the index n = s where / s (p) = 1 is the constant form factor. 



3 In particular, since the form factors do not depend on frequency, we assume that the frequency 
dependence of the vertex function is described well by the singular frequency lo- 
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The evolution of R is given by the particle-hole graphs, which unlike the particle- 
particle graphs remain bounded at all scales. In the beginning of the flow D is still 
small, so we neglect all particle-hole graphs proportional to D. We further neglect 
the remaining direct particle-hole graphs since they cancel excactly for a constant 
vertex function, hence remain small for a nearly constant vertex function. (The 
initial condition is V\ = U, that is, constant.) Then the flow equation for R, given 
by crossed particle-hole graphs only, is solved by R(q,l — q,q') = ip(q' — q) — U 
with ip(l) = U[l + U^Q,)]' 1 where $ ph (0 = J dp G(p)G(p + I) is the particle-hole 
bubble. 

Although <& ph (l) < (for Iq = 0), the denominator of if) remains nonzero for 
a sufficiently small initial coupling U since the particle-hole bubble is bounded in 
case of a curved and regular Fermi surface. So is a regular (smooth for an 
appropriate cut-off) and symmetric function that can be expanded in the form 
factors we are interested in. Again, the frequency dependence of ip is not taken into 
account here. Neglecting possible remainders of this expansion, the flow equation 
for the superconducting boson propagator reads 

Dmn{l) =<?n Yl [ D ™a(l) ~ *£'(0 [^WQ " K'n] , (20) 

a.a'eXsc 

where b aa , = f i>{q' - q) i ^ / a (q)/ a ' (q')- 

If we only take into account two form factors I sc = {s, 1} with /s(p) = 1 (s-wave) 
and /i(p) = cosp x — cosp y (d x 2_ y 2-wave), then b s i = b± s = 0, so b is diagonal as is 
<J?pp (0) in this case. Then the boson propagator D mn (0) is diagonal as well and the 
flow equation (I20p factorizes for I = such that the right hand side is always positive. 
The coefficient 6 SS is of order U and positive, so D ss is suppressed to zero. That is, 
s-wave superconductivity is always suppressed by the initial repulsive interaction. 
On the other hand &n is of order U 2 and negative, so -Dn(O) grows without bounds. 
Therefore, particle-hole fluctuations induce an attractive d-wave pairing interaction, 
which will be essential for the study of superconductivity in the Hubbard model 
at Van Hove filling. Here we have not included triplet superconductivity, which 
can become dominant as well. In general, attractive interactions in the Cooper 
channel are induced, so our method does capture the Kohn-Luttinger effect. The 
approximations made are valid provided the scale is not too low. 

In summary, the proposed decomposition of the interaction vertex is well under- 
stood for a curved and regular Fermi surface, where there is no Van Hove singularity 
and where Umklapp scattering is irrelevant. Namely, if the expansion in form factors 
is diagonal, we gave a clear argument which form factors have to be included in the 
expansion, such that the remainder of the expansion is negligible. In the general 
case diagonality holds approximately. The influence of subdominant particle-hole 
fluctuations can become important if there is no strong attractive interaction already 
initially. 

However, if Umklapp scattering and a Van Hove singularity in the density of 
states are present, particle-hole fluctuations can diverge by themselves and the mix- 
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ing between particle-particle and particle-hole channels is strong. This case is not 
analytically understood yet and treated numerically in the following. 



5 Q— Scheme and Integration over High Scales 

Although the renormalization group idea is independent of the type of regularization 
used to define the flow (provided the regularization satisfies some minimal conditions, 
e.g. that it makes the flow equation well-defined), its choice can become a delicate 
matter once approximations are used. In particular, we restrict the flow to the 
symmetric phase and apply the one-loop truncation. Parts of the vertex function 
will grow strongly at small scales, indicating an instability to an ordered (symmetry 
broken) state. Since we are using a weak coupling scheme, we have to stop the 
flow before all modes are integrated out, that is, the full model is not recovered 
in the flow. Therefore, the one-loop flow confined to the symmetric phase may 
depend on the type of regularization, and one has to compare the results of different 
regularization schemes. 

Since we are especially interested in the interplay between d-wave superconduc- 
tivity and ferromagnetism, we have to choose a regularization that does not suppress 



small momentum particle-hole fluctuations 13] . In the Hubbard model at Van Hove 
filling the density of states diverges logarithmically at low single-particle energies. 
At zero temperature and without scale regularization, the particle-hole bubble at 
zero momentum transfer is given by minus the density of states. Therefore, the 
particle-hole bubble at zero momentum transfer and zero temperature should di- 
verge logarithmically at low energy scales in the presence of a Van Hove singularity. 
However, for example, this is not the case for a Fermi surface cut-off which, at all 
nonzero scales, suppresses the small-momentum particle-hole excitations completely 



(see the discussion in 12| ) . For such a Fermi surface cut-off the particle-hole bubble 
at zero momentum transfer is zero for all finite scales at temperature zero. Even 
for nonzero temperature this means an artificial suppression of small momentum 
particle-hole processes down to the very lowest scales. Instead of a Fermi surface 
cut-off we multiply the bare propagator with the regulator function 

< 21 » 

which is independent of momentum but depends on frequency. The frequency f2 is 
the scale parameter we use to generate the RG flow and it replaces the A we used in 
the general discussion in Section [2j Since possible zeros of the denominator in the 
integrand of the bubble integration are canceled, the bubbles are regularized for all 
temperatures if the scale Q > 0. Small momentum particle-hole processes are not 
artificially suppressed to lowest scales with this regularization since the particle-hole 
bubble, now formed with the fi-dependent propagator G(p) = [ipo — e(p)] Xn(p), 
has the right asymptotic scaling behavior at temperature zero and Van Hove filling, 
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that is $ ph (0) ~ log f2. The particle-particle bubble $ pp (0) ~ (log f2) 2 diverges faster 
for the same parameters, as expected. 

In the limit O — > oo the regulator function (|2ip vanishes, hence the vertex 
function is the initial on-site interaction U of the Hubbard model. For f2 — > 0, 
Xn(p) — * 1) i- e - the regulator is removed. The regularization is mild, but it suffices 
to make all loop integrals in the flow equation converge. Thus it is an admissible 
scheme. For brevity we call this regularization the O-scheme. 

Because perturbation theory converges for large Q 2l|, [27|, we fix a scale £Iq 



and perform the integration over the degrees of freedom with propagator (1 — 
Xq (p)) [*Po — e (p)] by perturbation theory. The result provides the initial condi- 
tion for the flow at scale = Oo- To second order in U this gives for the Hubbard 
model in terms of a vertex function as in equation ([3]), 

Vn (pi,P2,P3) = U- U 2 $ pp ( Pl + P2 ) - U 2 % h ( P3 - Pl ) + O(U/Q ) 3 , (22) 

where the first term is the initial Hubbard repulsion. The second term is the particle- 
particle bubble, which generates an attractive s-wave superconducting channel. It 
will always be dominated by the repulsion U. The third term in (|22p comes from the 
crossed particle-hole graph. The contribution from the direct particle-hole graphs 
cancels out because the Hubbard interaction is on-site. 

We choose the initial condition for the boson propagators at scale f^o a s 

D ssfi (l) = U 2 <5> PP {1) (23) 
M SS)0 (0 = K aafi (l) = -U 2 % h {l) > 

for the constant s-wave form factor / s (p) = 1. The symmetric part of D SSi q(1) is 
positive as well, as is L> SSi o(^o = 0) !)• This assignment of the result of perturbation 
theory to the initial boson propagators is not unique. For example, attracting parts 
of the particle-hole term could be absorbed in a superconducting d-wave term. But 
the assignment (|23p is consistent with the definition of the three channels in Section 
[31 which do not mix in second order. Furthermore, the Kohn-Luttinger effect is very 
small in perturbation theory. The attractive coefficient of d x 2_ y 2 -superconductivity 
in an expansion of U 2 ^ pb (ps — p\) is bounded by ~ 0.01U 2 even for low scales and 
remains finite if the scale regularization is removed. Therefore the ambiguity of (|23p 
will not be important. In fact, we have numerically checked different distributions of 
the interaction terms in the initial conditions. We choose £lo > U big enough such 
that the results presented in Section [5] do not depend on the initial conditions for the 
boson propagators. The latter are negligible compared to interactions generated in 
the flow by the local repulsion term linear in U. So for large enough ° ur results 
are independent of the particular choice of Qq. 



6 RG Equations of the Boson Propagator Flow 

In order to derive the RG equations for the boson propagators we insert the ex- 
pansion of the three channels (|12p and (|14p into the RG equation §S§ and project 
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the right hand side of it according to the definitions of the boson propagators 
and (fT5"j) . Similarly, equations for the remainder terms R sc , ^m, and R K can be 
obtained. 

In case of a curved and regular Fermi surface we discussed in Section U how 
to separate singular from regular parts in the superconducting channel if the other 
channels are neglected. That gave a condition which form factors should be included 
in the expansion of the superconducting channel. Furthermore, as also analyzed in 
Section HJ we developed a good understanding of the integration of all one-loop 
graphs down to intermediate scales in this case. Therefore we argue that, with 
the right choice of form factors, the remainder terms can be dropped because the 
channels are described well even by only a few boson propagators. 

In presence of Van Hove singularities and relevant Umklapp scattering, however, 
the mixing of the channels is very hard to control. Here we rely on the results of 
[f| 0, [3) 2) 13] and choose a combination of form factors corresponding to the main 



instabilities found there. A detailed study of the remainder terms is left for future 
work. 

The flow equations derived in this section are valid for arbitrary index sets I sc 
and 2mk 5 assuming they both contain the index s that stands for the constant form 
factor /s(p) = 1- If the remainder terms are dropped, then the flow equation for 
the superconducting boson propagator reads for n, m e 2~ sc 

D mn {l) = <J m \ d/x(p, l-p)[ D rna {l)f a {\ - p) - U6 mS - a%(p, I)] (24) 

[ J2 D an (l)f a (t-p)-U6 nS -a s n c (p,l) , 

"■else 

where d//(p, I ± p) = dp-^G(p)G(l ±p) is the bubble integration. The initial local 
repulsion U will suppress s-wave superconductivity D ss . The contribution of the 
magnetic and forward scattering channels is given by 



X 



l) = \l 'rS^-q) E fi><*?)M-*F) 



6,6'6X MK 

(2 + a m )M w (p - q) - a m K w (p - q) 



qo=f 



For example, the magnetic boson propagator M ss induces an attracting interaction 
for d x 2_ y 2-wave superconductivity as described in Section HI The corresponding 
diagrams of the flow equation (|24p are plotted in Figure [1] for singlet symmetry 
c-m = +1 in a symbolic way. The graphs in the brackets are the superconducting, 
initial Hubbard repulsion, magnetic, and forward scattering interaction respectively. 
The square on the right hand side means that any two of the graphs in the bracket 
are put beside one another and connected to form one-loop diagrams. This generates 
three types of graphs, namely direct, vertex correction, and box graphs (see Figure 
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2). The latter two are not of the form of a boson propagator mediating between two 
fermion bilinear s. We expand them in a sum of terms of that form (given by our 
ansatz for the interaction) and drop the remainder term. To get the coefficients in 
the expansion (fl2l) . we apply the projections equation. 



5snwx=p sc 

w mn 



Y ,¥ 

+ 2 

A A 



Figure 1: The graphical representation of the flow in the singlet superconducting 
channel. 

In order to illustrate how the square in Figure [U is applied, we state two ex- 
amples. Taking two superconducting interactions and connecting them with two 
fermion propagators such that they form a loop, gives the direct graph in Figure 
[2|a). If the loop integration is performed, it is again of the form of two Cooper 
pairs interacting via a boson. By connecting a superconducting interaction with a 
magnetic interaction a vertex correction graph arises, given in Figure E^b). This 
is not of the demanded form, so the projection is nontrivial. A box diagram, for 
example, is formed by connecting two magnetic interactions, see Figure [2|c) . 

(a) 




V 



v v 



A. A ,A 



(b) (c) 

Figure 2: Three examples of arising graphs if the square in Figure [T] is taken: (a) 
the direct graph, (b) a vertex correction graph, and (c) a box graph. 

In the magnetic channel the flow is described by 



S) = ~ I dfi( P ,p + l)[ ]T M mb {l)f b {p + \) + U5 mS + aZ{p,l)] (25) 

[ Y, M bn(l)fb(p + k) + U6 nS + cC(p,l)], 
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where 




^/ m (q+^) ~2a m D (m ,(p-q)f a (£p + l)f a ,(Ep) 



a,a'£l sc 



+ E h{^)h>{\+^) M bb ,{p- q )-K w (p-q) 




Again the initial Hubbard interaction contributes only to the local s-wave coupling 
part. It will drive the boson propagator M ss . In contrast to that, the superconduct- 
ing s-wave boson propagator D ss will screen the local magnetic interaction M ss . 
Because this effect is quite substantial, especially if U is not very small, D ss is im- 
portant and must be included in the flow, even though it does not become singular. 
The influence of d x 2_ y 2 superconductivity with form factor /i(p) = cosp x — cosp y 
depends strongly on the boson momentum 1. It will suppress M ss (l) near 1 = 0, 
that is, suppress ferromagnetism, and enhance M ss (l) near 1 = tt = (jr, 7r), that is, 
antiferromagnetism, since /i(p + n) = — /i(p). The symbolic representation with 
graphs of the flow equation ([25]) is given in Figure O 



-,2 




Figure 3: The flow in the magnetic channel. 



Finally the flow equation for the forward scattering channel is given by 




^GXmk 



mS 



(26) 



X 



E K bn {l)fb(v + h)-U5 mS -al{p,l) 



with 




X /a'(l+ E i H ) 



+ E A( E ?)A'd + E | a ) 3M w (p-q)+K w (p-q) 




fc,f/eZMK 
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Similarly to s-wave superconductivity K ss is suppressed by U. In analogy to the 
Kohn-Luttinger effect a Pomeranchuk boson propagator K\\ is induced by magnetic 
correlations. Graphically the equation is sketched in Figure HI 




Figure 4: The flow in the forward scattering channel. 

We finish this section with some remarks on the just derived RG equations. 

1. Instead of the Hubbard model consider a mean-field model. So set U = and 
let the initial condition for the boson propagators be given by D mn fl{l) = D^n^ifi 
and similarly for M and K. Then the flow equations recover mean-field theory in 
the thermodynamic limit. Only the direct graphs contribute in the thermodynamic 
limit since the a's acquire an additional factor. 

2. Observing the flow equations, we find that the diagonal parts of the boson 
propagators at frequency zero are monotonically increasing. This is no contradic- 
tion to previous ./V-patch schemes where the measured susceptibilities did not need 
to be monotone. For example, the one-particle irreducible part of the connected 
expectation value 

-lj dpf m ^-p)Hq)^P(l-qy, j d P f n (±-p)iP(q)^(l-q)\ , 

which corresponds to the singlet superconducting susceptibility with f n (p) = f n (—p), 
is here given by D mn (l) plus a contribution from the other channels. If the latter 
prevail, the slope of the subdominant susceptibility can change sign even for n = m 
and I = 0. 

3. The flow equations stated so far are quite general. They are valid for arbitrary 
sets of form factors X sc and X M k- Also, no assumptions have been made concerning 
the frequency and momentum dependence of the boson propagators. For a numerical 
solution we specify this momentum dependence in the next section. For given 2"sc 
and T MK this leads to substantial simplifications. 

7 Momentum Dependence of the Boson Propagators 

While we have chosen particular functions for the form factors, the frequency and 
momentum dependence of the boson propagators is not specified yet. If this depen- 
dence is not constrained, the RG equation for the boson propagators still remains a 
system of integro differential equations which is hard to solve. For the following nu- 
merical calculation we neglect the frequency dependence of the boson propagators. 
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Together with the frequency independent form factors, this corresponds to a vertex 
function that does not depend on frequency at all. The right hand side of the flow 
equations is evaluated at zero boson frequency, which gives the main contribution. 

Furthermore we approximate the boson propagators by two-dimensional step 
functions. That is, their momentum dependence is discretized by dividing momen- 
tum space into segments on which the boson propagators are constant. Then the 
RG equation is equivalent to a system of x ordinary differential equations where x 
is the number of boson propagators times the number of segments. Note that the 
segments chosen here are neighborhoods of a single point because the boson prop- 
agators have point singularities. So they are different from the patches around the 
Fermi surface used in TV-patch studies. 

Thus the computational effort is reduced compared to Fermi surface iV-patch 
schemes, where the frequency dependence is also neglected and x ~ N 3 with N 
the number of patches for one fermion momentum. In many iV-patch studies N 
is reduced by projecting the momentum dependence onto the Fermi surface. This 
is not done here, we allow a general discrete momentum dependence of the boson 
propagators. It will turn out in the RG flow that the singular momentum dependence 
of the boson propagators is mainly determined by the fermion one-loop bubbles. 
Therefore, contrary to iV-patch schemes, we can choose the size and alignment of 
the segments to gain a more accurate approximation by step functions. 

For small scales and low temperature the bubbles can develop strong peaks at 
transfer momentum 1 = and 1 = tt = (tt, 7r) (or at 1 = fr — S for incommensurate 
Umklapp scattering). Possible singularities of the boson propagators are closely 
related to the corresponding direct graphs, which have the same momentum de- 
pendence as the bubbles. Therefore we separate the momentum dependence of the 
boson propagators into two parts 

B(l) = (\l x \ + \l y \ < vr) + - tt)1 (\l x - tt\ + \l y - tt| < tt) (27) 

for B = D, M, or K representing the different boson propagators and l x , l y £ [0, 2tt) 
and periodically continued elsewhere. We now approximate B^ a \l), where a = 
or 7T, by step functions for \l x \ + \l y \ < tt with l x ,l y € (— tt, tt]. Singularities will 
appear for small 1 = (l x ,l y ) only. However, since we do not choose !/ « 1, the 
magnetic local propagators Mgg\l) driven by U are not negligible for 1 away from 
zero. Likewise, the propagators (1) and (1) are of order U for all momenta. 
Therefore the step functions have to cover full momentum space and not just a small 
neighbourhood of 1 = 0. 

Symmetries of the boson propagators can be used to reduce the number of step 
functions necessary. If the boson propagators are diagonal, then they obey the 
symmetries B^ a \l x ,l y ) = B^ a \-l x ,-l y ) = B^(l y ,l x ) where a = or tt. In fact, 
for these symmetries to hold, the boson propagators are allowed to consist of block 
matrices for special form factors, such as /i and f^. Changing to polar variables 
lx = pcoscp and l y = psin</>, we divide the radial coordinate p € [0,^] into n 
intervals as can be seen in Figure [5j We choose smaller segments for small p, that is, 
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where a > 1 is a parameter 



for k = 1, ... n the segments are [( 71' (n)° 71 
for choosing the radial size of the smallest segment. The boson propagators are 
evaluated at the left border of each interval. Using the symmetries stated above, 
the angular variable <p needs only be discretized in [—§,§]■ This range is divided 
into m homogeneous intervals, where the boson propagators are evaluated at the 
midpoints. The k = 1 radial segment has no angular dependence imposed. Since 
we want to cover \l x \ + \l y \ < w, one extra segment is needed for the difference of 
the disk to the square. This more inaccurate approximation for large momenta is 
justified since the singularity of the bubbles develops at small momenta near or jr. 
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Figure 5: Schematic plot of the segment sizes for n = 4, m = 3, and a = 2. The 
circles mark the spots where the boson propagators are evaluated. The momentum 
dependence of the boson propagators on full momentum space is obtained by using 
symmetries. 



In the numerical calculation presented in the next section, we find that the 
magnetic and forward scattering boson propagators have sizeable contributions on a 
larger radial momentum range than the superconducting boson propagators, which 
have a strong peak for very small momenta only. This can also be observed in 
a comparison of the particle-particle and the particle-hole bubbles. Therefore we 
allow different a's for these two cases. 

Omitting the frequency dependence and parameterizing the momentum depen- 
dence with step functions simplifies the flow equations drastically for given form 
factors. The Matsubara frequency sum over the loop frequency po can then be 
calculated explicitly. Also, by trigonometric identities the q-integrations in the pro- 
jections in a^fi, a^, and become sums of terms, whose dependence on the loop 
momentum p can be stated analytically. Thus only the momentum loop integral 
has to be performed numerically to compute the flow. 
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8 Numerical Results at Van Hove Filling 



We have numerically solved the flow equations for the boson propagators with 
different sets of form factors. Already a flow with only a constant form factor 
/s(p) = 1 in all channels plus standard d :E 2_^2-superconductivity with form factor 
/i(p) = cos Px — cos p y in the superconducting channel gives reasonable results. We 
have checked that including the higher d x 2_ y 2 harmonic /s(p) = cos 3p x — cos3p y in 
the superconducting expansion produces only very small changes. Here we present 
results that are obtained by expanding all channels in the same set of form factors 
X sc = I MK = {a, 1}. Since 3> P p(a) = $pf(a) = for a = or it, the boson propaga- 
tors are diagonal at momentum 1 = or it. Furthermore they are diagonal for all 
\l x \ = \l y \ and this does not change in the flow. However, away from \l x \ = \l y \ off- 
diagonal terms appear, which at first are neglected here since the singular structure 
of the boson propagators lies on these lines. At the end of this section we present 
results where off-diagonal boson propagators are included. 

The density is set to Van Hove filling, that is, when varying the next to nearest 
neighbor hopping —t' the particle number is changed to fit /J, = At' . With this choice 
of parameters the Hubbard model at weak coupling has a strong tendency towards 
ferromagnetism. Since we want to compare our results with those of the temperature 



RG flow [12j|, where ferromagnetism was first found in a one-loop RG method for 
the Hubbard model, we choose U = 3 as wellfl Smaller couplings are harder to treat 
numerically. 

While in principle our method can treat all temperatures, for simplicity we here 
concentrate on temperature zero, which is not accessible within the temperature RG 
flow. However, starting at Qq = 15, we find that for all choices of t' , at least one 
of the boson propagators becomes almost singular at a point. In an approximation 
(p 2 + m^) -1 for the boson propagator, one would say that the mass of the 
boson tends to zero. However, the propagators are not well-approximated by this 
simple form, except maybe at very small momentum (depending on the scale £1), and 
this is the main reason why we use a numerical method to capture the momentum 
dependence. The flow to strong coupling observed in previous fermionic RG studies 
is related to this pointlike singularity, in that certain approximations correspond 
to setting p = 0, which hides the fact that the interaction only becomes strong at 
points. 

Although in our study, there is no flow to strong coupling in that sense, a true 
singularity in the interaction function introduces a significant change in power count- 
ing, which seems best to be captured by stopping the fermionic flow at a certain 
scale f2 c , performing a transformation to exchange bosons, and then attempting to 
continue to the symmetry-broken phase. In our present study, we only consider the 
flow to the scale fi c , which we call the critical scale, and leave the other steps to 
future work. We stop the flow before the maximum of any boson propagator reaches 
-Bmax = 20. The thus defined critical scale, which roughly corresponds to a critical 



4 As noted before, all quantities are measured in units of t = 1. 
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temperature in the temperature flow, is plotted in Figure over different hopping 
—t'. The flow to strong coupling of a boson propagator is interpreted as an insta- 
bility to a corresponding ordered state. However, this can only give a qualitative 
picture of the actual ground state phase diagram. For the latter the flow has to 
be continued in the symmetry broken phase. So for example, the specific choice of 
B max is arbitrary. However, it does not effect the qualitative features of Figure EJ 
For .B max = 16, for example, the critical scale is at most about 15% higher in the 
magnetic regions and even less in the superconducting region. Also, the range of —t' 
where superconductivity dominates shrinks by less than 1%. 
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Figure 6: The critical scale in dependence on hopping —t' at temperature zero and 
Van Hove filling for initial coupling U = 3. The instabilities of the Landau Fermi 
liquid are determined as antiferromagnetism (AF), d x 2_ y 2 superconductivity (dSC), 
and ferromagnetism (FM). The figure is further explained in the text. 

As described in Section [7] we approximate the momentum dependence of the 
boson propagators by step functions, which are constant on segments in momentum 
space. Even for huge segments with m = 1 and n = 4 qualitatively consistent results 
are obtained, although the crossover values of —t' between different instabilities as 
well as the critical scale are shifted compared to Figure El where we have used m = 5 
and n = 14. For small and big — t' the critical scale shows little dependence on dif- 
ferent segment sizes, especially on the angular resolution. However, in the crossover 
region between superconductivity and ferromagnetism for intermediate — t 1 the an- 
gular dependence becomes important, although the qualitative features of Figure 
remain unchanged. Plotting the momentum dependence obtained in the flow, it can 
be seen that only the boson propagators B^' away from small —t' have a signifi- 
cant angular dependence. Those boson propagators have four identical maxima that 
move away from 1 = as —t' is increased. In an angular sector that contains a max- 
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imum, they slowly increase on a high plateau for small increasing radial momentum. 
For higher radial momenta in the same sector they fall off rapidly after they passed 
the maximum, see also Figure [TIT a) at the end of this section. On the other hand, in 
between angular sectors that contain a maximum, there are sectors where the value 
of the boson propagator decreases from 1 = on with increasing radial momenta. 
Consider for example Figure [7(a), where the antiferromagnetic boson propagator 
Mgjp i s plotted at scale £l c an d parameter —t' = 0.3. At this hopping the sym- 
metric state is unstable towards d x i_ y i superconductivity, which is induced by the 
antiferromagnetic boson propagator. Neglecting the angular momentum dependence 
ofM s 7 results in a slightly higher J7 C for superconductivity. 
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Figure 7: Momentum dependence of the boson propagators (a) Mgg (antiferromag- 
netism) and (b) (d x 2_ y 2 superconductivity) at hopping —t' = 0.3 and critical 
scale O c . The dependence is obtained by an approximation with step functions 
characterized by m = 5, n = 14, and a = 2 in (a), and a = 3 in (b). Although the 
peak in (b) is very sharp, the small segments chosen at small momenta describe it 
by more than one constant. 



Generally, the magnetic and forward scattering boson propagators have sizeable 
contributions for a wide range of radial momenta, while the superconducting boson 
propagator D\\ has a very strong peak at 1 = 0, especially in the region where 
superconductivity is dominant. This can be seen in Figure[7Jb) where D\\ is plotted 
at O c and hopping — t' = 0.3. In order to adapt to these different momentum scales 
we have chosen a sc = 3 for all superconducting boson propagators and a MK = 2 for 
all magnetic and forward scattering boson propagators. 

Now we discuss Figure [6) For small — t' the boson propag ator M!£\l) at 1 = 
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grows strongest, indicating an instability towards an antiferromagnetic phase. Since 
the flow is stopped at a relatively high scale O c ~ 0.1, perfect nesting-like effects 
are not restricted to —t' = 0. If —t' is increased between 0.08 < —t' < 0.23 the 
antiferromagnetic boson propagator is still the leading term. However its maximum 
is away from 1 = 0, that is, it has four identical maxima, similar to Figure[7{a). This 
reflects a tendency to incommensurate antiferromagnetic order, indicated by open 
squares in Figure El 

For —t' € [0.23, 0.34] the leading term is the superconducting boson propagator 
(1) at 1 = 0, which represents d x 2_ y 2~wave superconductivity. As discussed in 
Section^ it is induced by antiferromagnetic correlations. In Figure EKa) the flow of 
D^(0) is plotted for —t' = 0.3 in dependence on the scale. First incommensurate 
antiferromagnetic correlations dominate and (0) is enhanced by them but stays 
small down to low scales. When a significant superconducting coupling is reached, 
the direct graph ~ (0) 2 <j>pp(0) contributes strongly such that D^(0) increases 
very fast and becomes the dominant coupling. The crossover between the antiferro- 
magnetic and the d x 2_ y 2 superconducting region is not sharp in the scheme applied 
here. Around —t' = 0.23 both tendencies grow very strongly, making it difficult to 
judge numerically which is the dominant instability. The overlap of the two regions 
in Figure [9^a) corresponds to the saddle point region found earlier in [15]. 

As —t' is increased the antiferromagnetic correlation, indicated by |$ ph (-7r)|, de- 
creases, while ferromagnetism is enhanced because |<& ph (0)| becomes bigger. Both 
affect the critical scale in the superconducting region, which is continuously drop- 
ping, especially when —t' gets close to the ferromagnetic region. This indicates that 
d x 2_ y 2 superconductivity is not only induced by antiferromagnetic correlations but 
also suppressed by ferromagnetic tendencies. 

For high —t' > 0.34 the leading instability is ferromagnetism, represented by 
the dominant flow of M^g (1 = 0). The critical scale rises again as —t' increases. 
However, the critical scale is still about one order of magnitude lower than the 
Stoner criterion suggests. The Stoner criterion for ferromagnetism is obtained by 
neglecting the particle-particle channel and given by £/|<£ ph (0)| = 1. Therefore 
superconducting fluctuations suppress ferromagnetism as well. However, outside 
the d-wave superconductivity region the superconducting d-wave boson propagator 
(1) remains small, see Figure E^a) for its maximum at 1 = 0. A significant 
contribution to this suppression comes from s-wave superconductivity D$g and 
so-called screening, compare Figure [TU1 

In summary, ferromagnetism and d-wave superconductivity are competing insta- 
bilities of the (t, £')-Hubbard model for relatively large —t' at Van Hove filling. In the 
calculation presented here the crossover between superconductivity and ferromag- 
netism takes place at —t' = 0.34 which is in good agreement with —t' = 0.33 found 
with the temperature RG flow [lil ] and the two-particle self-consistent Monte Carlo 
approach [§]. This is also roughly the value of — t' where ^ P h(0) = <I> ph (7r). However, 
compared to the temperature RG flow, the suppression between these competing 
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Figure 8: The dominant flows for two different values of —t' in the superconducting 
region. The peak of the boson propagator at momentum zero evolves slowly 
but reaches the value 20 first. According to the definition of the critical scale Q c , 
the flow is stopped and all boson propagators are manually put constant for scales 
below J7 C . In (a) —t' = 0.3 is well inside the superconducting region while in (b) 
—t' = 0.34 is the crossover to the ferromagnetic region. Here both instabilities grow 
strongly. Their mutual suppression results in a lower critical scale. In both figures 
antiferromagnetic fluctuations are dominant for intermediate scales, inducing d x 2_ y 2 
superconductivity. 



tendencies is weaker in the approximation used here. Although for low and high —t' 
the critical scale £l c and the critical temperature from 12j is roughly the same, the 
temperature RG flow finds critical scales that are two orders of magnitude lower in 
the crossover region between superconductivity and ferromagnetism. Especially, the 
results obtained here do not suggest a quantum critical point separating these two 
instabilities in the phase diagram. Instead, in the crossover region both processes 
grow strongly like in the crossover between antiferromagnetism and superconduc- 
tivity. However, the crossover is much sharper as can be seen in Figure El^a). The 
superconducting boson propagator Dy± remains small as soon as —t' enters the fer- 
romagnetic region, while in the antiferromagnetic region superconducting tendencies 
persist, such that this crossover region is broader. 

At present, it is not clear to us where this difference to [13] comes from, as 
both schemes contain rather different approximations. The RG equations for the 
boson propagators imply that their diagonal parts are monotonically increasing in 
the flow. Already a small inexactness, such as the omission of the remainder terms 
or the discretization of the momentum dependence, can prevent an exact cancella- 
tion on the right hand side of the flow equation for and M^. On the other 
hand, the projection to the Fermi surface, used in [12] but not here, is expected to 
enhance superconducting tendencies. This could possibly lead to a stronger sup- 
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pression. Furthermore, the fi-scheme used here allows a clearer definition of the 
initial conditions than the temperature scheme. It would be interesting to perform 
an iV-patch analysis with the proposed Q-regularization ([2Tj) . 




Figure 9: Overview of the maximal value of the boson propagators at the critical 
scale, which is plotted as a thin dotted line. In (a) the three dominant instabilities 
are plotted, that is antiferromagnetism, d x 2_ y 2 superconductivity, and ferromag- 
netism. The maxima of some sub leading boson propagators are plotted in (b). The 
kinks in the plot indicate the crossovers to another dominant instability. For large 
—t 1 the discretization scheme for the momentum dependence becomes insufficient 
for the antiferromagnetic boson propagator, whose four maxima move towards the 
edge of its support. Due to the large segment sizes there, maxiMgg^(l) jumps dis- 
continuously if the maxima leave the support and another segment contains the new 
maxima, compare Figure [71 



So far we have not discussed all boson propagators computed in the flow. The 
s-wave forward scattering boson propagators and and the s-wave super- 
conducting boson propagators Ds$ and D^s cannot become singular since they are 
dominated by the initial local repulsion U. In Figure Mjo) their final values are 
plotted over hopping —t'. They are biggest in the crossover region between super- 
conductivity and ferromagnetism, indicating the amount of screening and the low 
critical scale. The boson propagator K-j^ , representing a possible Pomeranchuk 
instability, is induced by magnetic correlations in the flow. However, it is not 
dominant in the parameter region considered here. The other boson propagators 
D^p , m[® , Mfi , and K-j^ remain smaller than 0.05 and hence are irrelevant here. 
Therefore, in order to further reduce computing cost, they could be omitted. In 
fact, only the boson propagators , Mq , and K$ do not satisfy the symmetry 
B( a \l x ,l y ) = B^ a \l x , —ly). If they are not computed in the flow, the angular seg- 
ments only need to cover the interval [0, This would halve the number of step 
functions necessary for the discretization of the momentum dependence described 
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in Section [71 

We have already mentoined that the effective on-site interaction changes al- 
though the initial on-site repulsion U of the Hubbard model is kept constant in the 
parametrization of the vertex function ([9]). Projecting to the on-site part gives the 
scale-dependent effective on-site coupling 

U ctt = U+f if^H \-DS (1) + I A4 a) (1) - (1)" , 

J\l*\+\ly\<* ( 27r ) ~ # L 2 2 J 

which is plotted in Figure [TU1 over hopping —t'. Although the critical scale is lowest 
in between the superconductor and the ferromagnet, the effective on-site coupling 
has not increased much in comparison to higher scales. This indicates that the 
amount of screening, especially in the crossover region between superconductivity 
and ferromagnetism, is substantial. 




Figure 10: The effective on-site coupling U cH is plotted in dependence of — t' for 
different scales in (a). The curves to different values of Q are drawn only over the 
interval where that is still bigger than the critical scale f2 c . In (b) U eS is evaluated 
at the critical scale, which is plotted as a dotted line for comparison. The kinks in 
(b) indicate the crossovers of distinct instabilities. 

In a second numerical computation we study the influence of offdiagonal boson 
propagators. In the superconducting channel we consider the boson propagators 
Z?ss\ £>ss\ £>ii\ and D$ . The first two are important to capture screening effects, 
the third describes d x 2_ y 2 superconductivity and the last is an offdiagonal boson 
propagator that couples s- and d-wave superconductivity. In the magnetic and 
forward scattering channel we only take into account the local s-wave boson propa- 
gators Mgg\ Mgg , , and K^. Compared to the previous numerical calculation, 
apart from the irrelevant boson propag ators D[f, M$ , M-ff, and K$ , only the 
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boson propagator , which indicates a possible Pomeranchuk instability, is omit- 
ted. However, remained relatively small in the first computation, see Figure [9j 
Note that for this choice of boson propagators there are no more offdiagonal terms 
than the one included. 

The flow equations imply that the diagonal boson propagators obey the addi- 
tional symmetry B^ a \l x , l y ) = B^ a \l x , —l y ). The offdiagonal boson propagator 
obeys this symmetry as well and is also symmetric under 1 — ► —1. However, it is 
antisymmetric under exchange of l x <-* L. Therefore, compared to the choice of step 
functions described in Section [71 we only need to consider segments with an angle 
4> = arctan j- in the interval [0, Since we still use the parameter values m = 5 
and n = 14, the angular resolution is higher in this second numerical computation. 
Again we computed the RG flow of the boson propagators for coupling U = 3 at 
temperature zero and Van Hove filling for various —t' € (0, |). Since no major 
process was omitted compared to the first numerical computation, no significant 
change could be observed in the figures presented. Furthermore, the maximum of 
the offdiagonal boson propagator remains smaller than 0.3 for —t' > 0.37 and 
smaller than 0.03 for —t' < 0.37. The offdiagonal boson propagator that couples s— 
and d-w&ve superconductivity can therefore be neglected in good approximation. 

The only difference in the two different numerical calculations is found in the 
shape of the antiferromagnetic boson propagator Mgf^ for relatively high —t'. We 
already commented on the insufficient description of large radial momenta of the an- 
tiferromagnetic boson propagator Mgg \ which is plotted in Figure [TT1 at —t' = 0.3 
in dependence on p in different angular sectors, in (a) for the first and in (b) for the 
second numerical computation of this section. Note that the angular degeneracy in 
(a) is due to the same number of angular segments in the larger angular interval 
[— f, f ]■ The major difference between both figures is that in (b) the boson propaga- 
tor has its maximal value at smaller radial momentum p. This can play a role in the 
search for an analytic parametrization of the antiferromagnetic boson propagator, 
although the plateau in the angular sector <p = for intermediate p seems to be its 
most important feature. However, note that the angle <fi = is not a center of an 
angular sector in (b). The maximal value of the boson propagator is found in the 
nearest angular segment, which contains <j) = but is centered at <f> = This could 
also be a partial explanation of the observed difference. 



9 Conclusion 

We have presented a novel parametrization of the one-loop one-particle irreducible 
RG equation for the four-point function and applied it to the two-dimensional Hub- 
bard model using a novel RG regularization scheme. The parametrization is based 
on the idea to separate dominant interactions from irrelevant remainders, where 
the latter do not influence the leading instabilities of the flow. This idea is imple- 
mented by writing the effective two-fermion interaction as a sum of dominant terms 
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(a) (b) 




Figure 11: The antiferromagnetic boson propagator M ss plotted over radial mo- 
menta for different angular sectors centered at d> = arctan y- . For the first nu- 
merical calculation shown in (a) the centers of the angular sectors take the values 
4>i = — j + (i — ^)jq and for the computation with an offdiagonal boson propagator 
in (b) the angular sectors are centered at <j>j = (j — 2)^0 • 



where fermion bilinears interact via exchange bosons. Each such term consists of 
a boson propagator and boson-two-fermion vertex functions. From the singular 
momentum structure of the RG equation we identify three different channels in the 
interaction. In these channels the boson-fermion vertex function is expanded in 
scale-independent form factors. We argued that only a few terms are needed in the 
expansions to describe the qualitative structure of the flow. The tails of the expan- 
sions are the above mentioned remainder terms, which are neglected. Although the 
scale dependent coefficient functions of the expansions are called boson propagators, 
they are kept as fermionic interaction terms. Their flow is computed by suitable 
projections of the fermionic RG equation. 

The application of this method to the (t, t')-Hubbard model at Van Hove fill- 
ing and temperature zero shows indeed that the essential qualitative structure of 
the one-loop RG is preserved. We are able to reproduce the leading weak cou- 
pling instabilities of the model, as they were found in previous RG studies using 
iV-patch schemes. In particular, our results are in good qualitative agreement with 
the temperature RG flow . We use a new frequency regularization, which like 
the temperature regularization does not artificially suppress ferromagnetism, but 
which allows to state the initial conditions of the RG flow more clearly. For an ini- 
tial repulsive on-site interaction U = 3 we find three distinct regions, characterized 
by the leading weak coupling instability and depending on next to nearest neighbor 
hopping, see Figure [U For small next to nearest neighbor hopping —t' antiferromag- 
netism dominates. As was found in the A-patch schemes, antiferromagnetic correla- 
tions induce a d^^-wave superconductivity instability, which becomes dominant 
for intermediate —t'. For high —t' the leading instability is ferromagnetism. In 
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between the d x 2_ y 2-wave superconducting and the ferromagnetic region the critical 
scale drops by two orders of magnitude. We argue that this is due to the mutual 
suppression of these opposing correlations. Compared to 12j], however, we find that 
this suppression is weaker here. In particular, our result does not suggest a quantum 
critical point between the two regions. Because we have not been able to trace this 
discrepancy to certain scattering processes, we regard this matter as open and leave 
it for further study. 

The proposed method is in general not restricted to temperature zero and Van 
Hove filling. Finite temperature only affects the calculation of the one-loop fermion 
bubbles, which then becomes more involved. Away from Van Hove filling triplet 
superconductivity plays a role for high —t' [12]. The flow equations stated in Section 
[6] can account for that. However, it is not clear that a single form factor suffices 
to describe dominant triplet superconductivity. The results of 12] suggest that at 
least three different triplet form factors would have to be taken into account. 

The benefit of the proposed parametrization of the vertex function is a reduction 
of the complexity of the full one-loop flow to some dominant terms. We believe 
that separating leading from subleading processes will help to gain further insight 
into the structure of the one-loop RG. However, we have not performed a detailed 
analysis of the remainder. The comparison with previous A-patch studies suggests, 
nevertheless, that we did capture the most important processes. A rigorous proof is 
particularly complicated at Van Hove filling due to the strong mixing of the particle- 
particle and the particle-hole channels. We gave clear arguments for determining 
the dominant terms in case of a regular and curved Fermi surface. 

Practically, this reduction of complexity results in lower computing cost com- 
pared to previous A-patch schemes. In the numerical implementation the momen- 
tum dependence of the boson propagators is discretized using step functions. The 
choice of segments, where the step functions are constant, can be guided by the form 
of the one-loop bubbles. This allows a more precise discretization than the general 
patching of the vertex function in A-patch schemes. In particular, no momentum 
dependence is projected onto the Fermi surface. Furthermore, the proposed decom- 
position gives possibly way to further improvement. If the momentum dependence 
of the one-loop bubbles can be parameterized in an analytical form, then it should 
be possible to extract a functional parametrization of the boson propagators from 
the flow equations, at least for small momenta. Deviations for large momenta, away 
from the maxima of the boson propagators, could be subject to another negligible re- 
mainder. However, the momentum dependence of the bubbles is difficult to describe 
since a naive power expansion is not sufficient. Nevertheless, while in the numerical 
calculation we have neglected the frequency dependence of the boson propagators, 
a similar procedure for small frequencies could take into account at least part of the 
frequency dependence of the vertex function. 

The proposed decomposition of the effective two-fermion interaction is of a 
form that suggests the decoupling of the fermion bilinears via multiple Hubbard 
Stratonovich transformations. The ambiguity of introducing boson fields is not 
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completely removed but at least reduced. This allows to continue the RG flow into 
the symmetry broken phase in a (partially) bosonized form. 

We acknowledge financial support from DFG grant SA 1362/1-2 and DFG re- 
search group FOR 723. 
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